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摘 要 : 基于 1956 一 2021 年 芍 勒 河上 游 昌 马 河 出 山口 昌 马 保水 文 站 的 实测 洪水 数据 及 流域 内 或 周 
边 鱼 儿 红 雨量 站 和 托 勒 气象 站 的 气温 和 降水 数据 ,采用 Sen'sslope 估 计 法 `Mann-Kendall 检验 法 ， 
抽样 分 析 了 年 最 大 洪水 (Annual maximum peak discharge,AMPD ) 系 列 、 超 定量 国 值 洪水 (Peaks- 
over-threshold ,POT3M ) 系列 和 基于 广义 极 值 (Generalized extreme value,GEV ) 分 布 不 同 级 别 洪 水 系 


2 


7 


的 变化 特征 。 结 果 表 明 :(1) 1956—2021 F R LB BB EK KH LK A 


为 显著 性 增加 的 趋势 ,其 中 AMPD 和 POT3M 系列 的 幅度 变化 趋势 通过 了 显著 性 水 平 为 0.01 的 检 


验 ,POT3M 系 列 的 频次 变化 趋势 通过 了 显著 性 水 了 


为 0.05 的 检验 ,基于 GEV 分 布 不 同 级 别 洪水 以 


小 洪水 为 主 , 且 总 洪水 频次 通过 了 显著 性 水 平 为 0.01 的 检验 。(2) 1987 年 以 后 气温 和 降水 呈现 波动 


性 增加 趋势 , 强 降 雨 和 升温 引起 的 冰川 .高山 积 雪 


和 频次 显著 性 增加 的 主要 原因 。 
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文章 编号 : 1000 - 6060(2023)08 - 1260 - 09( 1260 ~ 1268) 


河流 洪水 是 人 们 最 关心 的 自然 灾害 问题 之 
一 。 根 据 联合 国有 关 部 门 的 统计 数据 ,1995 一 2015 
年 洪水 是 全 世界 所 有 记录 自然 灾害 中 最 频繁 的 事 
件 , 占 所 有 灾害 的 43% ,全 球 每 年 平均 损失 估计 为 
1040x10° USD 。 洪 水 产生 的 影响 和 不 同 区 域内 洪 
水 的 变化 特征 已 经 引起 了 广泛 的 关注 ,并 成 为 世界 
范围 水 文学 领域 许多 研究 的 焦点 2 。 如 洪水 的 幅 
度 和 频次 的 增加 或 减少 与 区 域 范 围 内 的 温度 上 升 、 
强 降 水 和 水 文 循 环 加速 等 水 文 要 素 有 关于 。 

从 已 有 研究 来 看 ,与 气候 变化 有 关 的 洪水 变化 
特征 已 经 成 为 许多 研究 的 热点 问题 ,主要 涵盖 洪水 
事件 的 量 级 .频次 和 时 间 变 化 的 影响 。 最 近 的 
研究 表明 ,欧洲 范围 的 气候 变化 不 仅 改变 了 洪水 的 
发 生 时 间 ,而 且 由 于 秋季 和 冬季 增加 的 降水 导致 了 
西北 欧 地 区 洪水 的 增加 ,相反 由 于 降水 减少 温度 
升 高 和 蒸发 量 增加 导致 东南 部 洪水 的 减少 ,已 
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冻 土 融 水 增加 等 原因 是 3 种 抽样 洪水 系列 幅度 


有 学 者 在 中 国 西北 的 干旱 和 半 干 旱地 区 也 开展 了 
一 些 研 究 , 如 Zhang 等 "发 现 20 世 纪 90 年 代 后 塔 里 
木 河 流域 的 洪水 规模 呈现 增加 的 趋势 ,特别 是 位 于 
高 纬度 的 河流 。 蒋 军 新 等 ”和 毛 炜 峰 等 ”发 现 天 
山 山 脉 的 典型 流域 ,如 阿克苏 河 、 托 什 干 河 、 库 玛 拉 
克 河 ` 玛 纳 斯 河和 乌鲁木齐 河 ,在 过 去 50 a, 随 着 气 
候 变 化 的 影响 ,洪水 量 级 明显 增加 。Wang 等 也 发 
现 1970 一 2019 年 发 源 于 祁连山 的 12 条 河流 洪水 以 
夏季 小 洪水 为 主 , 洪 水 频次 呈 上 升 趋势 ,东部 和 西 
部 地 区 存在 差异 ,东部 地 区 在 减少 ,西部 地 区 在 增 
加 。 芷 勒 河 正好 位 于 祁连山 和 河西 走廊 的 西部 , 研 
究 表明 上 游 出 山口 昌 马 堡 水文 站 多 年 平均 径流 量 
约 10.3x10: mx ,上 游 的 来 水 量 对 极度 缺 水 的 中 下 
Vet PM KARAS RRB KES, mI tbh 
ERE AS A kU A Us A SY ee "ELL 
区 高 大 山系 分 布 较 多 ,冰川 和 积 雪 发 育 ,径流 主要 
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由 降水 和 冰雪 融 水 形成 “。 因 此 ,在 气候 变化 背景 
下 ,分 析 和 探究 玖 勒 河 历 史 时 期 的 洪水 变化 规律 ， 
对 中 下 游人 民 财 产 安全 和 经 济 社会 发 展 具有 重要 
意义 。 

人 研究 基于 过 去 60 多 年 下 勒 河上 游 昌 马 河 出 山 
口 昌 马 堡 水 文 站 的 长 序列 实测 洪水 数据 ,采用 年 最 
大 洪水 (Annual maximum peak discharge, AMPD ) 系 
列 、 超 定量 病 值 洪水 (Peaks-over-threshold,POT3M) 
系列 和 基于 广义 极 值 (Generalized extreme value， 
GEV) 分 布 不 同 级 别 洪水 系列 ,分 析 了 1956 一 2021 
年 昌 马 堡 水 文 站 洪水 系列 的 幅度 和 频次 变化 趋势 ， 
并 试图 探究 洪水 的 变化 特征 与 气候 变化 的 归 因 关 
系 。 以 期 为 当地 防汛 抢险 水利 工程 建设 及 调度 管 
理 . 水 资源 管理 等 提供 决策 依据 。 


1 研究 区 概况 


研究 区 概况 如 图 1 所 示 。 玖 勒 河 发 源 于 祁 连 
山脉 的 岗 格 尔 首 合力 岭 , 河 源 高 程 4737 m, TRE 
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长 583 km, , 集 水 面积 4.13x10' km ,地 理 位 置 介 于 
92°54'~99°14'E ,38°36'~41°34' N 之 间 。 从 河源 到 昌 
hie tH LW (B SRC) Aaa al EDE, RZ 
昌 马 河 ,河流 长 347 km, 集 水 面积 约 为 1.10x10* km’, 
地 理 位 置 介 于 96°42’ ~99°00'E,, 38°18’ ~39°54'N 之 
间 。 干流 从 河源 区 开始 由 东南 向 西北 穿行 于 托 
来 南山 与 下 勒 南山 之 间 ,两 山 分 水 岭 高 程 一 般 在 
4000m 以 上 ,其 中 臣 勒 南山 的 最 高 峰 是 室 召 吉 
WE ,高 程 为 5808 m, 是 祁连山 脉 的 最 高 峰 。 上 游 海 
拔高 程 介 于 2080~5808 m ,平均 海拔 3944 m, HHI E 
主要 由 大 雪山 oie BS LU FER FA LL A A aT e 
组 成 ,两 侧 山 区 地 势 高 峻 、 地 形 陡峭 ,谷地 地 形 相对 
平缓 " 沾 。 下 勒 河上 游 由 于 海拔 高 ,气候 寒冷 ,成 为 
现代 山岳 冰川 和 冻 土 发 育 区 。 冰 川 主要 分 布 于 政 
勒 河谷 两 侧 的 大 雪山 a FA LL PFE OR Fa LL PE 
4248 m 以 上 的 区 域 ,覆盖 面积 达 469.5 km’, 24 r5 Yii 
域 面积 的 4.3%"“*。 人 研究 区 总 体 呈 现 高 原 山 地 气候 
特征 。 


(b) Bia TT pits 


图 1 研究 区 观测 站 点 分 布 


Fig. 1 Distribution of observation stations in the study area 


2 数据 与 方法 


2.1 数据 来 源 

选取 足 勒 河 上 游 昌 马 河 出 山口 昌 马 保水 文 站 
为 洪水 观测 站 点 ,流域 内 鱼 儿 红 雨 量 站 和 周边 托 勒 
气象 站 为 气象 观测 站 点 。 鱼 儿 红 雨量 站 位 于 上 游 
区 间 有 连续 降水 观测 资料 ,本 文 参考 河源 区 周边 托 


勒 气象 站 的 数据 对 比分 析 洪 水 变化 特征 对 降水 和 
气温 变化 的 响应 。 水 文 站 和 气象 站 信息 如 图 1 和 
表 1 所 示 。 洪 水 数据 采用 昌 马 保水 文 站 1956 一 
2021 年 洪水 观测 值 。 昌 马 保 水 文 站 和 人鱼 儿 红 雨量 
站 的 数据 来 源 于 《中 华人 民 共 和 国 水 文 年 鉴 》, 托 勒 
气象 站 的 数据 来 源 于 中 国 气 象 局 信息 中 心 (http:// 


data.cma.cn )。 
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表 1 研究 区 观测 站 点 基本 信息 


Tab.1 Basic information of observation stations in the study area 


河流 观测 站 点 观测 项 流域 面积 /10’ km? 海拔 /m 经 度 纬度 数据 系列 
HLA THY 昌 马 堡 水 文 站 流量 1.0961 2080 96°51" 39°49" 1956—2021 年 
鱼 儿 红 雨 量 站 降水 - 3030 97°06' 39°26' 1956—2021 年 
FCM TA 气温 、 降 水 - 3367 98°25! 38°48’ 1956—2021 年 


2.2 研究 方法 

2.2.1 Sen’s slope 估 计 RH Sen’ s slope Hit at 
析 玻 勒 河 昌 马 堡 水 文 站 AMPD 洪水 系列 POT3M 潜 
水 系列 和 基于 GEV 分 布 不 同 级 别 洪 水 系列 的 幅度 
和 频次 变化 趋势 ,以 及 鱼 儿 红 雨 量 站 和 托 勒 气象 站 
的 气温 和 降水 变化 趋势 。 


Yi 
p= median =A 了 > 天 (1) 


式 中 :6 为 正 值 或 负 值 表示 上 升 或 下 降 趋势 ,其 值 的 
大 小 表示 趋势 的 斜率 倾斜 程度 ;和 x 为 时 间 序 列 
数据 , 且 x 和 代表 时 间 j7 和 的 数据 值 。 
2.2.2 Mann-Kendall 检验 法 采用 Mann-Kendall 
检验 法 (M-K 检 验 ) 对 玻 勒 河 昌 马 保水 文 站 AMPD 潜 
水 系列 .POT3M 洪水 系列 和 基于 GEV 分 布 不 同 级 别 
洪水 系列 的 幅度 和 频次 变化 趋势 ,以 及 鱼 儿 红 雨量 
站 和 托 勒 气象 站 的 气温 和 降水 变化 趋势 进行 分 析 。 
M-K 检 验 法 是 一 种 广泛 应 用 于 水 文 和 和 气象 要 素 
时 间 序 列 趋势 分 析 的 非 参 数 检验 方法 。M-K 检 验 法 
所 选 的 样本 不 必 先 假设 服从 正 态 分 布 ,也 不 受 个 别 奇 
异 特征 数 值 的 影响 ,其 详细 表达 式 见 参考 文献 1。 
2.2.3 POT3M 类 水 系列 抽样 方法 POT3M 洪 水 系 
列 中 ,洪水 阅 值 的 选择 目前 还 尚 无 公认 的 客观 方 
法 '。 部 分 学 者 通过 试 算 确 定 平均 每 年 发 生 洪水 
次 数 介 于 2.4~3.0 之 间 对 应 的 洪水 流量 为 阔 值 ,本 研 
究 最 终 通过 试 算 确定 洪水 发 生 的 频次 为 2.7 场 .a '。 
为 了 使 抽样 洪水 洪峰 的 独立 性 假设 可 以 成 立 ,采用 
了 Lang 等 提出 的 洪水 独立 性 判别 方法 ,同时 考虑 
了 在 中 国 西 北 干 量 和 半 干 旱地 区 应 用 的 POT 抽样 
方法 "中 。 在 本 研究 中 ,所 有 抽样 的 洪水 洪峰 在 持 
续 时 间 (D) 和 中 间 最 小 流量 (Qi,) 方 面 都 符合 以 下 2 
个 条 件 : 


. >5+log(A) a) 
2 


Ons Fmin(Q,,Q.) 
式 中 : 刀 为 两 场 连 续 洪 水 洪峰 流量 之 间 的 洪水 持续 
时 间 (d) ;4 为 集 水 面积 (km) ;0 为 两 场 连续 洪水 


洪峰 流量 0, 和 0Q, 之 间 的 最 小 流量 (ms ')。 

2.2.4 基于 GEV 分 布 不 同 级 别 洪水 系列 抽样 方法 FE 
于 GEV 分 布 拟 合 AMPD 洪水 系列 频率 分 布 曲 线 , 采 
用 K-S 检 验 法 和 PPCC 检验 法 进行 拟 合 优 度 检验 和 
相似 程度 评价 ””。 检 验 通过 后 计算 洪水 理论 重 现 
期 并 确定 其 相应 的 洪水 流量 浆 值 ,根据 不 同 重 现 期 
的 洪水 流量 国 值 将 洪水 等 级 分 类 为 小 洪水 .中 洪水 
和 大 洪水 。 对 于 干旱 和 半 干 旱地 区 的 玲 勒 河 来 说 ， 
每 年 的 洪水 数量 较 少 ,参考 《水 文 情报 预报 规范 
(GB/T 22482-2008)》 中 对 洪水 等 级 的 划分 ,确定 相应 
洪水 等 级 分 为 :小 洪水 ( 重 现 期 <5 a) ;中 洪水 (5 a< 
重 现 期 <20 a) ;大 洪水 ( 重 现 期 >20 a)。 

GEV 分 布 函数 为 : 


etre re io (ho 


(3) 
式 中 : x 为 洪水 序列 数据 ; é 为 形状 参数 ; o 为 尺度 
参数 ; u 为 位 置 参 数 。 当 #<0 时 ,服从 Weibull 分 布 ; 
当 E=0 时 ,服从 Cumbel 分 布 ; 当 E>0 时 ,服从 Frechet 
分 布 。 本 研究 采用 极 大 似 然 估计 法 对 GEV 分 布 进 
行 参数 估计 ,其 似 然 函数 为 : 


1 


L(é,o,1) =-n Ino - $j eN i [ + x 


i 


Sf 4 (4) 


式 中 :n 为 时 间 序 列 长 度 ;x; 为 时 间 i 的 洪水 序列 数据 。 
3 结果 与 分 析 


3.1 AMPD 洪水 系列 变化 特征 

1956—2021 年 琉 勒 河 昌 马 保水 文 站 AMPD 洪 水 
系列 变化 特征 如 图 2 和 表 2 所 示 。 年 最 大 洪峰 流量 变 
异 系数 值 为 0.43, 极 值 比 为 8.05, 年 际 变 化 幅度 极 大 ， 
最 大 洪峰 流量 出 现在 1999 年 ,为 718 ms ,而 最 小 洪 
峰 流 量 出 现在 1956 年 ,为 89.2 m-s MREZE 
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表 2 1956—2021 年 琉 勒 河 昌 马 堡 水 文 站 AMPD 洪水 系列 变化 特征 
Tab.2 Characteristics of AMPD flood series at Changmabao hydrological station in the Shule River from 1956 to 2021 


平均 洪峰 最 大 洪峰 


最 小 洪峰 


河流 “水 文 站 极 值 比 
a 流量 /ms 变异 系数 “流量 ax.s HME SLEE BEms 出 现年 份 ”与 多 年 均值 比 
me Baye 322 0.43 718 1999 2.23 89.2 1956 0.28 8.05 
120 


距 平 百分率 /% 


y=1.066x—2119, R°=0.219, P<0.01 


150 
1950 1960 1970 1980 1990 2000 2010 2020 2030 
年 份 


TE: AMPD 为 年 最 大 值 抽样 方法 的 洪水 系列 。 
图 2 1956 一 2021 年 踊 勒 河 昌 马 堡 水 文 站 AMPD 洪水 系列 
局 度 变 化 趋势 
Fig.2 Trend in the magnitude of AMPD flood series at 


Changmabao hydrological station in the Shule River from 
1956 to 2021 


势 和 5 a 平 滑 曲线 来 看 ,AMPD 洪水 系列 为 增加 趋 
势 ,线性 倾向 率 为 1.066,1956 一 1987 年 为 保持 平稳 
期 , 昌 大 多 数 年 份 AMPD 值 小 于 距 平 值 ,1988 一 2021 
年 为 增加 趋势 , 晶 1999 年 以 后 大 多 数 年 份 AMPD 值 
明显 大 于 距 平 值 。M-K 检 验 中 AMPD 洪 水 系列 通过 
了 显著 性 水 平 0.01 的 检验 ,为 显著 性 增加 趋势 。 
3.2 POT3M 洪水 系列 变化 特征 
3.2.1 POT3M 类 水 系列 幅度 变化 POT3M 洪 水 系 
列 幅度 变化 趋势 如 图 3 所 示 。 从 线性 变化 趋势 和 5 a 
平 请 曲线 来 看 ,洪水 幅度 为 增加 趋势 。 与 AMPD Hk 
水 系列 相似 ,1956 一 1990 年 为 保持 平稳 期 ,有 小 幅 
波动 , 除 1981 年 为 正 距 平 外 ,其 余 大 部 分 为 负 距 平 ; 
1991—2021 年 为 波动 增加 趋势 。M-K 检 验 中 POT3M 
潜水 系列 的 幅度 变化 通过 了 显著 性 水 平 0.01 的 检 
验 , 为 显著 性 增加 趋势 。 
3.2.2 POT3M 洪水 系列 频次 变化 

POT3M 洪 水 系列 频次 变化 趋势 如 图 4 所 示 。 从 
线性 变化 趋势 和 5 a 平滑 曲线 来 看 ,洪水 频次 为 增加 
趋势 。1956 一 1985 年 有 小 幅 波 动 , 大体 为 平稳 期 ; 
1986—2021 年 为 波动 性 增加 趋势 。 与 POT3M 洪 水 
系列 的 幅度 变化 不 同 , M-K 检 验 中 频次 变化 通过 了 
0.05 的 显著 性 检验 ,为 显著 性 增加 趋势 。 


| Ml POT3M 一 一 5 a 平滑 曲线 ------ 线性 拟 合 


EFAN% 


y=0.669x—1334, R°=0.211, P<0.01 


20 
1950 1960 1970 1980 1990 2000 2010 2020 2030 
年 份 


注 :POT3M 为 超 定量 阔 值 抽样 方法 的 洪水 系列 。 下 同 。 
图 3 1956—2021 年 巩 勒 河 昌 马 堡 水 文 站 POT3M 洪 水 系列 
幅度 变化 趋势 
Fig.3 Trend in the magnitude of POT3M flood series at 


Changmabao hydrological station in the Shule River from 
1956 to 2021 


Ml POT3M 一 一 5 a 平滑 曲线 --…--- 线性 拟 合 


12 y=0.064x-124, R°=0.157, P<0.05 
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图 4 1956—2021 年 巩 勒 河 昌 马 堡 水 文 站 POT3M 洪水 系列 
频次 变化 趋势 
Fig.4 Trend in the frequency of POT3M flood series at 
Changmabao hydrological station in the Shule River from 
1956 to 2021 


3.3 基于 GEV 分 布 的 洪水 变化 特征 

3.3.1 不 同 级 别 洪水 的 重 现 期 冰 值 计算 ”采用 GEV 
分 布 拟 合 玖 勒 河 昌 马 保 水 文 站 AMPD 洪水 系列 , 拟 
合 参 数 和 拟 合 优 度 检验 值 见 表 3。 从 表 3 的 结果 分 
析 可 知 ,形状 参数 (&) 大 于 0, 服 从 Frechet 分 布 。K- 
S 检 验 值 为 0.109, 小 于 显著 性 水 平 0.05 的 临界 值 
0.196, 拟 合 精 度 良好 。PPCC 检验 值 为 0.95 ,表明 
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昌 马 保水 文 站 AMPD 洪水 系列 服从 GEV 分布 的 
Frechet 分 布 , 且 高 度 相 关 。GEV 分 布 对 分 析 计 算 
玻 勒 河 昌 马 堡 水文 站 的 重 现 期 闪 值 有 很 好 的 可 靠 
性 。 同 时 采用 MATLAB 软件 程序 计算 1956 一 2021 
年 跑 勒 河 昌 马 保水 文 站 的 重 现 期 闽 值 ,结果 如 表 4 
所 示 。 


表 3 1956—2021 年 琉 勒 河 昌 马 堡 水 文 站 AMPD 
系列 拟 合 参 数 与 检验 值 
Tab.3 Fitted parameters and test values of the AMPD 
series at Changmabao hydrological station in the 
Shule River from 1956 to 2021 


水 文 站 形状 参数 ”尺度 参数 ”位 置 参 数 ” 拟 合 K-S PPCC 

(£) (o) (u) 分 布 ”检验 检验 

昌 马 堡 0.018 107 258 Frechet 0.109 0.95 
分 布 


表 4 1956—2021 FRAI SS 7k ASMA 
Tab. 4 Threshold values of the return period at 
Changmabao hydrological station in the Shule River 
from 1956 to 2021 


重 现 期 赣 值 /as” 


河流 水 文 站 
la Sa 20a 
Ia 昌 马 堡 98 420 585 


15 fr (a) 小 洪水 
mame 小 洪水 — 5 a 平滑 曲线 ------ 线性 拟 合 
y=0.069x—132, R=0.294, P<0.01 
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SEAS 
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5 场次 /1987 一 2021 年 


频次 / 场 
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3.3.2 不 同 级 别 洪水 频次 变化 ， 采 用 GEV 分 布 分 
别 计算 大 洪水 ( 重 现 期 >20 a) .中 洪水 (5 a< 重 现 期 
<20 a) 和 小 洪水 ( 重 现 期 <5 a) 的 频次 变化 趋势 (图 
5) 。1956 一 2021 年 洪水 发 生 的 频次 主要 以 小 洪水 
JE , Sen’ s slope 估 计 倾 向 率 为 0.085, H. M-K 检 验 
中 总 洪水 的 频次 通过 了 0.01 的 显著 性 检验 ,洪水 频 
次 变化 为 显著 性 增加 趋势 。 
3.3.3 不 同年 代 际 的 洪水 频次 变化 ”从 图 6 可 以 看 
出 ,1956 一 2021 年 琉 勒 河 昌 马 保水 文 站 不 同年 代 际 
洪水 频次 总 体 为 增加 趋势 。 各 年 代 相 比较 而 言 ， 
1950s—1960s 洪水 发 生 频 次 较 少 ,1970s 开始 增加 ， 
且 1970s 一 2000s 为 保持 平稳 期 ,2010s 达到 洪水 发 
生 频 次 最 大 值 。 总 体 来 看 , 玻 勒 河 洪水 发 生 频 次 演 
变 趋势 为 由 少 向 多 转换 。 
3.4 降水 和 气温 变化 特征 

由 图 7 和 图 8 可 知 ,1956 一 2021 年 鱼 儿 红 雨量 
站 和 托 勒 气象 站 的 年 降水 量 以 及 托 勒 气象 站 的 气 
温 总 体 呈 现 显 著 性 增加 的 趋势 。 鱼 儿 红 雨 量 站 和 
托 勒 气象 站 的 多 年 平均 降水 量 分 别 为 140.3 mm 和 
302.4 mm, Sen’ s slope 估 计 倾 向 率 为 6.9 mm- (10a)! 
和 15.3 mm = (10a) '; 托 勒 气象 站 的 多 年 平均 气温 
为 -2.4% ,Sen’s slope 估 计 倾 向 率 为 0.35 °C (10a) '。 


6r (b) 中 洪水 
5 场次 /1956 一 1986 年 
4 16 场 次 /1987 一 2021 年 
X 
x 
EN 
0 
1950 1960 1970 1980 1990 2000 2010 2020 2030 
年 份 
15- ® 总 洪水 


大 洪水 — 5 aP HHR ------ 线性 拟 合 
中 洪水 

127 小 洪水 

y=0.085x—164, R°=0.324, P<0.01 


0 
1950 1960 1970 1980 1990 2000 2010 2020 2030 
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图 5 1956—2021 年 芍 勒 河 昌 马 保水 文 站 不 同 级 别 洪水 频次 变化 趋势 
Fig.5 Trend in the frequency of different levels of floods at Changmabao hydrological station in the Shule River from 1956 to 2021 
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图 6 1956 一 2021 年 踊 勒 河 昌 马 保水 文 站 
不 同年 代 际 洪水 频次 
Fig.6 Frequency of interdecadal floods at Changmabao 
hydrological station in the Shule River from 1956 to 2021 


(a) 托 勒 气象 站 
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(b) 鱼 儿 红 雨 量 站 


y=0.689x—-1371, R=0.091, P<0.05 
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图 7 1956 一 2021 年 托 勒 气象 站 和 鱼 儿 红 雨量 站 
年 降水 量变 化 趋势 
Fig.7 Trends of annual precipitation at Tuole 
meteorological station and Yu’ erhong rainfall 
station from 1956 to 2021 


由 图 9 可 知 , 托 勒 气象 站 单 场次 降水 量 >10 mm 的 降 
水 日 数 和 年 平均 降水 强度 都 为 增加 趋势 。 由 图 7~9 
分 析 可 知 ,1987 一 2021 年 降水 量 和 气温 为 正 距 平 ， 
日 为 持续 增加 态势 。 
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图 8 1956—2021 年 托 勒 气象 站 气温 变化 趋势 
Fig. 8 Trend of temperature change at Tuole meteorological 
station from 1956 to 2021 


(a) 降水 日 数 


故国 降水 量 >10 mm ”一 一 5 a 平滑 曲线 ------ 线性 拟 合 
y=0.055x-101, R=0.123, P<0.01 
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(b) 降水 强度 
6 [i 降水 强度 一 5 a 平滑 曲线 -…-- 线性 拟 合 


y=0.028x—-53.7, R?=0.656, P<0.01 


降水 强度 /mm*d- 


0 
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图 9 1956—2021 年 托 勒 气象 站 降水 要 素 变化 趋势 
Fig.9 Trends in precipitation elements at Tuole 


meteorological station from 1956 to 2021 


4 讨论 


TT EDF KIB es Fe LL D ye E, 
从 河源 区 至 昌 马 堡 水 文 站 区 间 ,海拔 高 差 达 2000 多 
米 , 区 间 流 域 地 形 .坡度 和 下 垫 面 等 都 有 差异 ,上 游 
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区 域 主要 受 高 原 季 风 和 西风 带 大 气 环 流 系统 的 交 
汇 影 响 ,属于 西风 带 气候 区 。 通 过 分 析 流 域内 鱼 儿 
红 雨 量 站 和 周边 托 勒 气象 站 的 降水 变化 可 知 (图 7、 
图 9) ,过 去 几 十 年 , 玻 勒 河流 域 的 降水 量 呈 持续 增 
加 趋势 ,尤其 是 夏季 降水 量 占 总 降水 量 比重 越 来 越 
KO) 。 已 有 研究 成 果 也 证 实 了 降水 量 的 变化 趋 
势 , 如 施 雅 风 等 2 时 在 21 世 纪 初 期 研究 发 现 ,从 
1987 年 开始 整个 西北 地 区 的 气候 巾 暖 干 向 暖 湿 转 
型 ,表现 为 温度 上 升降 水 量 增 加 、 冰 川 消融 量 和 径 
流量 连续 多 年 增加 ; 蓝 永超 等 ”人 徐 浩 动 等 "研究 
发 现 , 玻 勒 河 上 游 降 水 量 总 体 为 增加 趋势 。 从 托 勒 
气象 站 的 气温 变化 趋势 可 知 (图 8),1987 一 2021 年 
气温 为 显著 性 增加 趋势 。 已 有 研究 发 现 ,1956 一 
2010 年 琉 勒 河 河 源 区 的 年 平均 气温 以 0.40 C- 
(10a) 的 倾向 率 显著 性 升 高 ”; 中 国 西北 干旱 区 气 
温 的 低温 极 值 在 1986 年 左右 突变 后 呈现 为 显著 性 
增加 趋势 ”在 温 升 背景 下 ,与 冰川 、 积 雪 和 冻 土 
等 冰冻 圈 要 素 相 关 的 下 垫 面 变化 ,对 玻 勒 河 等 西部 
寒 区 流域 产 流 机 制 已 产生 了 显著 的 影响 ”i。 第 一 
次 冰川 编 日 (1956 一 1983 年 ) 和 第 二 次 冰川 编目 
《2005 一 2010 年 ) 数 据 的 相关 研究 成 果 守 ”显示 p 
勒 河 冰川 覆盖 面积 和 冰 储 量 为 660 条 河 509.87 km? 
左右 。 随 着 温度 逐渐 升 高 ,1960 年 以 来 琉 勒 河 冰 川 
融 水 呈现 增加 趋势 邱 。Liu 等 外人 研 究 发 现 , 在 
1960s 一 2010s 期 间 气 候 持 续 变 暖 影响 下 ,冰川 融 水 
对 玻 勒 河 的 融 水 占 比 达到 42.2%。 对 比分 析 鱼 儿 红 
雨量 站 和 托 勒 气象 站 的 降水 和 气温 要 素 与 昌 马 堡 
水 文 站 洪水 发 生 的 幅度 和 频次 ,结果 为 正 相 关 关 
系 ,变化 趋势 大 体 一 致 。 


5 结论 


本 文 基于 1956 一 2021 年 政 勒 河上 游 昌 马 河 出 
山口 昌 马 保水 文 站 的 实测 洪水 数据 和 相关 气象 数 
据 ,分 析 了 AMPD 洪水 系列 .POT3M 洪 水 系列 和 基 
于 GEV 分 布 不 同 级 别 洪水 系列 的 变化 趋势 及 影响 
因素 。 主 要 结论 如 下 : 

(1) 3 种 抽样 洪水 系列 中 AMPD 系列 和 POT3M 
系列 的 幅度 和 频次 为 显著 性 增加 趋势 ,基于 CEV 分 
布 不 同 级 别 洪水 系列 以 小 洪水 为 主 , 且 总 洪水 频次 
为 显著 性 增加 趋势 。 

(2) 鱼 儿 红 雨 量 站 和 托 勒 气象 站 的 数据 表明 ， 


1987 以 后 气温 和 降水 呈 显著 性 波动 增加 趋势 , 强 降 
十 和 升温 引起 的 冰川 高山 积 雪 和 冻 土 融 水 增加 等 
原因 是 3 种 抽样 洪水 系列 幅度 和 频次 显著 性 增加 主 
要 驱动 因素 。 

琉 勤 河上 游 山区 的 洪水 变化 特征 受气 候 变 化 、 
冰冻 圈 要 素 变化 和 下 垫 面 等 原因 的 影响 十 分 复杂 ， 
本 文 仅 通过 数理 统计 方法 进行 了 简单 的 趋势 分 析 ， 
其 洪水 模拟 和 影响 机 制 有 竺 后 期 进一步 深入 人 研究 。 
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Characteristics of flood changes in mountainous areas of 
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Abstract: Based on measured flood data for 1956—2021 from the Changmabao hydrological station at the outlet 
of the Changma River in the upper reaches of the Shule River, northwest China, and on temperature and precipita- 
tion data from Yu’ erhong rainfall station and Tuole meteorological station in the vicinity of the basin, the varia- 
tion characteristics of annual maximum peak discharge (AMPD), peaks over threshold (POT3M), and floods 
based on the generalized extreme value (GEV) distribution at different levels were sampled and analyzed using 
Sen’ s slope estimator and the Mann-Kendall test. The results showed the following: (1) Both the magnitude and 
frequency of floods at the Changmabao hydrological station during the studied period significantly increased, 
with AMPD and POT3M magnitudes showing highly significant increases (P<0.01) and POT3M frequency show- 
ing significant increases (P<0.05). Based on the GEV distribution, floods were mainly small floods and the fre- 
quency of total floods was highly significant (P<0.01). (2) Temperature and precipitation showed a fluctuating in- 
creasing trend after 1987, and heavy precipitation and increased meltwater from glaciers, alpine snow, and perma- 
frost due to warming were the main causes of the observed significant increases in flood magnitude and frequen- 
cy in the three sampled series. 
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